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We study theoretically the superfluidity and stability of a Bose-Einstein condensate (BEC) whose 
interatomic scattering length is periodically modulated with optical Feshbach resonance. Our nu- 
merical study finds that the properties of this periodic BEC are strongly influenced by the modula- 
tion strength. When the modulation strength is small, only the Bloch waves close to the Brillouin 
zone edge suffer both Landau and dynamical instabilities. When the modulation strength is strong 
enough, all Bloch waves become dynamically unstable. In other words, the periodic BEC loses its 
superfluidity completely. 

PACS numbers: 03.75.Kk, 05.30.Jp, 67.10.Ba 



I. INTRODUCTION 

Even though superfluidity is one of the most impor- 
tant macroscopic quantum phenomena, it could only be 
found in liquid helium before 1995. Since the realiza- 
tion of Bose-Einstein condensation in atomic gases in 
1995, we have now in experiment another superfluid, the 
Bose-Einstein condensate (BEC) [1]. This new super- 
fluid shares many interesting properties with superfluid 
helium, such as critical velocity [2] and quantized vor- 
tices [3|. At the same time, there are also some inter- 
esting properties unique to this gaseous superfluid; for 
example, there is no roton excitation in BECs. In partic- 
ular, since a gas is easily compressible, one can modulate 
its density with an optical lattice to create a periodic su- 
perfluid. The properties of this periodic superfluid have 
been studied extensively both theoretically ^] and exper- 
imentally [il . Dynamical instability, which is absent in a 
homogenous superfluid, was discovered and found to play 
a dominant role in destroying superfluidity in a periodic 
superfluid 0, |5[. 

Now another type of periodic superfluid can be cre- 
ated: with optical Feshbach resonance (OFR) [ej, the in- 
teratomic interaction (or scattering length) of a BEC can 
be modulated periodically in space with laser beams [t^. 
This was already demonstrated in experiment [8]. This 
BEC with periodically modulated interaction (PMI) is 
different from the usual BEC in an optical lattice, where 
the atoms feel an external periodic potential (PP). The 
most important difference between these two periodic 
BEC systems is that the BEC with PMI has no lin- 
ear periodic counterpart. As a result, the widely used 
single-band approximation for a periodic system appears 
not applicable for a BEC with PMI. Furthermore, while 
the BEC in PP has a Mott-insulator phase d, [io|, the 



BEC with PMI should not have the Mott phase. In this 
work we study the superfluidity and other related physi- 
cal properties in a BEC with PMI. 

We focus on the case where the laser beam for OFR 
is applied only along one direction. In this case, the 
BEC can be described by the following Gross-Pitaevskii 
equation 
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\/H + [Vi + V2 cos{2kzx)]\^f^ , (1) 
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where Vi and V2 are positive parameters that can be 
tuned experimentally by changing laser power and de- 
tuning p-lsl; m is the atomic mass and kz is the wave 
number of the laser beam. There have already been some 
theoretic^ attempts to find soliton solutions in such a 
system [ijj, [l2| . In this work we study the superfluidity 
of this periodic system. This is equivalent to examining 
the stability of a flow described by a Bloch wave [l3|. 
The Bloch wave solutions of this BEC system for the 
lowest Bloch band are found numerically. Their Landau 
and dynamical stabilities are examined by computing the 
Bogoliubov excitations. We find that when the periodic 
modulation strength V2 is small, only Bloch waves close 
to the Brillouin zone edge are unstable. When the mod- 
ulation strength becomes large enough, all Bloch waves 
in the lowest band become dynamically unstable. This 
means that the periodic modulation of the scattering 
length can cause a BEC to lose its superfluidity com- 
pletely. When this happens, the periodic Bose system 
is neither a superfluid nor a Mott insulator; the BEC 
may collapse into many solitons as suggested by an early 
study [13. In contrast, for a BEC in PP, the Bloch state 
near the Brillouin zone center is always stable no matter 
how strong the periodic modulation is [4]. 

The paper is organized as follows. In Sec. II, we present 
the basic theoretical framework within which the BEC 
system is treated. In Sec. Ill, the lowest two Bloch bands 
for this periodic BEC system are presented and their 
physical meaning is discussed. In Sec. IV, we study the 
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Landau instability and dynamical instability of the Bloch 
waves and discuss the superfluidity of this system. In Sec. 
V, the results in the previous section are discussed in the 
context of the experiment. The paper is summarized in 
Sec.VI. 



II. THEORETICAL FRAMEWORK 



the system energy. This energetic instability is related to 
the superfluidity of the system and we call it the Landau 
instability as this is the essence behind Landau's the- 
ory of superfluidity. The Landau instability is sometimes 
called a thermodynamical instability. 

We also consider the dynamical evolution of the system 
after the perturbation Scpk^q. For such a small perturba- 
tion, the dynamical Eq.(j2j) can be linearized and becomes 



We consider the one-dimensional case, where the scat- 
tering length is modulated only in one direction by op- 
tical Feshbach resonance and the lateral motion of the 
BEC can be either ignored or confined. In this case the 
Gross-Pitaevskii equation becomes one-dimensional 



1 d'^ 

+ (ci + C2 cosx)|V^|^V^ . (2) 



Here the energy unit is 4:h'^k'^/m, and the length unit is 
l/2kz' The wave function ip is in units of ^/no^ where 
no is the averaged BEC density; Ci = and C2 = 

^^2^2 ' For convenience, we call ci the uniform strength 
and C2 the modulation strength. 

We are interested in the superfluidity and stability of 
a flow in this system. For a homogeneous BEC, the flow 
is described by a plane wave. For a periodic system, the 
flow is represented by a Bloch wave. The Bloch wave has 
the form = e*^^(/9/c(x), where (^k{x) is of period 27r, 
it satisfies the time-independent Gross-Pitaevskii equa- 
tion 
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+ (Cl + C2 C0SX)|V^PV^ . 



(3) 



We use the numerical method proposed in Ref. [15| to find 
these Bloch waves. 

Once a Bloch wave solution is found, its stability is 
examined. We add a small perturbation to the Bloch 
wave ^k{x)^ 



S(pk,q{x) = Uk{x, ^)e'^^ + vl{x, q)e 



-iqx 



(4) 



where q is in the range [—1/2,1/2] and represents the 
mode of perturbation. The energy deviation caused by 
the perturbation is 



dx{ul,vl)Mk{q) 



Uk 
Vk 



(5) 



where 

Mu{q) 
with 

m = 



C{k + q) {ci + C2 cos x)ipl , , . 



- jj.^ 2{ci^C2 COS x)\(pkf . (7) 



If the matrix Mk{q) has negative eigenvalues, it means 
that there are some perturbations S(pk,q{x) that can lower 



dt V Vk 



Vk 



I 
-/ 



(8) 



If matrix aMk{q) has complex eigenvalues, the system 
becomes unstable and will damp during dynamical evo- 
lution. We call this instability dynamical instability. 

Landau instability and dynamical instability can be 
discussed in a more coherent way within the framework 
of Bogoliubov excitations. The eigenvalues of matrix 
(jMk{q) can be divided into two groups, phonon modes 
and antiphonon modes; only phonon modes are physical 
and are usually called Bogoliubov excitations [15| . It can 
be proved that matrix Mk{q) having a negative eigen- 
value is equivalent to aMk{q) having phonon modes of 
negative energy [16]. Therefore, the system has Landau 
instability when some of the phonon modes have nega- 
tive energies and it has dynamical instability when some 
of the phonon modes are complex. As the phonon modes 
are related to the superfluidity of a Bose system [l^ , it is 
clear from this perspective that Landau instability and 
dynamical instability are clearly related to each other 
and are just two ways of destroying superfluidity. As we 
show, both instabilities are present in this periodic BEC. 



HI. NONLINEAR BLOCH BANDS 

Because of the periodic modulation of interatomic in- 
teraction, this BEC has Bloch wave solutions. These 
Bloch waves can in general be found by numerically solv- 
ing Eq.([3|). However, near the edge of the Brillouin 
zone(/c ^ 1/2), we can use a two- mode approximation 
and assume the Bloch state is of the form (pj^{x) ^ 
ae'Ck-h)^ + ^e^Ck+h)^ (|ap + |6p = 1, |^| <C 1 and 
k = k — 1/2). Plugging the trial wave function into 
Eq.(|3|), we can get a quartic equation [l8| 

+ 2gw^ + {g^ - - l)w^ - 2gw - g'^ = , (9) 

where g = Ci/c2, h = k/c2^ w = l/(2a6), and the chem- 
ical potential is /i = 1.5ci + 0.5c2{w -\- 1/w). It can be 
shown that near /c = 0, Eq.(|9]) has only two real solu- 
tions when ^ < 1 (ci < C2); it has four real solutions 
when g > 1 {ci > C2). This means that the Bloch band 
of this nonlinear periodic system has a loop structure at 
the edge of the Brillouin zone when g > 1. This is con- 
firmed by our numerical computation, as shown in Fig(TJ 
In Fig(T](b), where ci > C2, the chemical potential ja has 
a clear loop at the edge of the Brillouin zone. 
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FIG. 1: The lowest two Bloch bands of a BEC with its in- 
teratomic interaction periodically modulated, C2 = 0.1: (a) 
ci = 0.05, (b) ci = 0.2. /i is in units of 4h^k%/m and k is in 
units of 2kz, where kz is the wave vector the laser beam. 



The loop structure in the Bloch band in Fig(T] (b) is a 
manifestation of superfluidity in the system. We can con- 
sider the homogeneous case C2 = 0. In this case, the BEC 
is a superfluid with critical velocity Vc = Now we 

slowly turn on the periodic modulation by increasing C2 
to a small value. This small periodic modulation can be 
regarded as a perturbation when ci is large. The Bragg 
scattering caused by this periodic perturbation should 
not destroy a superflow moving with velocity v = 1/2, 
which is the velocity at the Brillouin zone edge when 
C2 ~ 0. The nonzero slope of the Bloch band at the 
edge is an indication of this robustness of superfluidity. 
When the periodic modulation becomes very strong as in 
Figdja), the Bragg scattering can eventually disrupt the 
superflow as indicated by the zero slope at the zone edge 
in Figdja). A similar loop structure has been found in 
many different systems [19]. Note that the loop struc- 
ture of the energy band has its interesting many-body 
counterpart, a net of narrow avoided crossings in energy 
levels iioi]. 



IV. SUPERFLUIDITY AND INSTABILITY 

For all the Bloch waves found in the lowest band, we 
have examined their superfluidity and stability by nu- 
merically computing their Bogoliubov excitations with 
(7Mk{q). The results are shown in the stability phase di- 
agram of FiglSl For a point (/c, q) in the figure, there are 
three possibilities: (i) If it is in the white region, it means 
that the Bloch wave is a local energy minimum rela- 
tive to the perturbation mode q. (ii) If it falls into the 
black area, the Bloch wave (/^/c is dynamically unstable 
relative to the perturbation mode q. Any perturbation 
containing mode q will cause the system to evolve dynam- 
ically away from state (pk with an exponential growth, 
(iii) If the point (/c, q) lies in the gray region, the Bloch 
wave (fk is not a local energy minimum but dynamically 
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FIG. 2: Stability phase diagrams of BEC Bloch states for 
different values of the uniform strength ci and modulation 
strength C2. A; is the wave number of the Bloch state and q is 
the wave number of the perturbation mode. In the white area, 
the Bloch state is a local energy minimum; in the shaded (light 
or dark) area, the Bloch state has negative excitation energy 
and Landau instability; in the dark shaded area, the Bloch 
state is dynamically unstable, k and q are in units of 2kz, 
where kz is the wave vector of the laser beam. 



stable relative to the perturbation mode q. Note that we 
have only plotted the results for /c > in Fig. [2] as the 
system is symmetric with respect to time reversal and 
the results are the same for k and —k. 

The stability phase diagram in Figl2fa) is very similar 
to the one for a BEC in PP [4]. However, as C2 increases 
or ci decreases, the phase diagram begins to have new 
features. The most prominent is that the black area (dy- 
namical instability) spreads into the region with k < 1/4 
and eventually reaches k = 0. For a BEC in PP, the 
black area is restricted in the region with A: > 1/4. It 
is also clear from the figure that the boundaries of the 
gray area and the black area coincide once the black area 
reaches /c = 0. This means that once the Bloch state at 
k = becomes unstable, it has both Landau instability 
and dynamical instability. In this case, the ground state 
of this system is no longer a Bloch wave. In contrast, the 
ground state of a BEC in PP is always a Bloch wave. 

There is a qualitative way to understand why the Bloch 
state at /c = becomes dynamically unstable. When 
ci < C2, the interaction between atoms becomes attrac- 
tive for some parts of the system. As C2 increases, a 
larger portion of the system becomes attractively inter- 
acting. Eventually, the overall interaction of the sys- 
tem, indicated by the averaged interaction over one pe- 
riod V = Jq^{ci +C2 cos x)\ilj\'^dx^ becomes negative. We 
expect that this is the underlying reason that the Bloch 
state at /c = becomes dynamically unstable. This is 
confirmed by our numerically computation. In FiglSl we 
have marked out the stability regions for the Bloch state 
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at /c = in the parameter space of ci and C2. We find 
numerically that the Bloch state at /c = is dynamically 
unstable in the black area in FiglH The solid(red) line 
is the dividing line between averaged positive interaction 
and averaged negative interaction. From Figl3]we can see 
the solid(red) line is lower than the stability boundary 
for the Bloch state. It means that the Bloch state is al- 
ways unstable when the averaged interaction is negative. 
However, this is not the only reason that the Bloch state 
becomes unstable as there is an area where the averaged 
interaction is positive while the Bloch state is unstable. 
This is yet to be fully understood. 

We now approximate the periodic BEC with overall 
negative interaction with a homogeneous BEC with at- 
tractive interaction strength c. For a homogeneous BEC, 
the flow is described by plane waves e*^^. In this case, 
the matrix crMk (g) is a 2 x 2 matrix. The phonon excita- 
tions are easily obtained. For the state with /c = 0, they 




FIG. 4: Stabilities of the Bloch waves in the lowest Bloch 
bands shown in Fig. 1. The Bloch states in the black area 
are dynamically unstable, and the Bloch states in the shaded 
area have Landau instability, C2 = 0.1: (a) ci = 0.05, (b) 
ci = 0.2. II is in units of 4h^k%/m and k is in unit of 2/cz, 
where kz is the wave vector of the laser beam. 



(10) 



Since c is negative, the phonon excitation is imaginary for 
small q. This is in fact the feature seen in Figsl2fc)-l2ff): 
the Bloch state at A: = is dynamically unstable against 
the perturbation of mode q = 0. This approximation re- 
sult further confirms that the averaged interaction plays 
a dominant role in this system of a BEC with PMI and is 
the key factor in stability and superfluidity of the BEC 
system. 




V. EXPERIMENTAL PERSPECTIVE 

The superfluidity and dynamical instability has been 
explored experimentally for a BEC in PP [5]. A similar 
experimental scheme can be used to study superfluidity 



FIG. 5: (Color online)The critical values of /cl and ko as a 
function of ci and C2. /cl is the critical Bloch wave number for 
the Landau instability and kn is the critical Bloch wave num- 
ber for the dynamical instability. The (red) lines are where 
kL and ko turn to zero. The two (red) lines in (a) and (b) 
are identical. kL and ko are in units of 2k where kz is the 
wave vector of the laser beam. 
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FIG. 3: (Color online) The stability phases of the Bloch state 
at /c = in the space of the interaction parameters ci and 
C2. The Bloch wave state is stable in the white area; it is 
dynamically unstable in the black areas. On the solid (red) 
line, the averaged interatomic interaction is zero; below it, 
the averaged interaction becomes negative; and above it, the 
averaged interaction is positive. 



and dynamical instability of this periodic BEC system. 
There exists no fundamental technical barrier. 

In a real experiment, there is no way to control the 
perturbation mode. The controlled perturbation and the 
uncontrollable noises in the experiment should contain 
all the possible modes of q. Therefore, for an experimen- 
talist, a Bloch state cpk is unstable if it is unstable against 
any of the perturbation modes. For the Bloch waves in 
the lowest band in FigdJ we have marked out their sta- 
bilities with gray and black shadings in FigHJ It is clear 
from the figure that the Bloch states near the Brillouin 
zone center are more stable and more Bloch states be- 
come stable as the uniform strength ci increases. 

There exist two critical values of k in FiglH kL and 
ko' A BEC Bloch wave cpk with k in the range [/cl, 1/2] 
has Landau instability; it has dynamical instability if its 
k is in the range [/ci:),l/2]. Both critical values kL and 
ko vary with ci and C2, and this dependence is shown in 
Fig IS As seen in FigEfa), kL increases with the increase 
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of ci and the decrease of C2. In FiglSJb), kn '^s discontin- 
uous at the edge of C2 = 0. This is due to that at C2 = 0, 
the system loses the periodicity and it has no dynamical 
instability. On the two (red) lines FiglH both /c^ and 
turn zero. The two (red) lines are the same. This reflects 
a fact that we already mentioned: when the Bloch wave 
at /c = has Landau instability, it is also dynamically 
unstable. 



VI. SUMMARY 

We have studied a BEG with its scattering length pe- 
riodically modulated in one direction. In this periodic 
BEG system, the flows are represented by Bloch waves 
and the energy has band structure. When the uniform re- 
pulsive interaction is larger than the periodic modulation 
strength, the band structure has a loop at the edge of the 
first Brillouin zone. We have also studied the stabilities 
of these Bloch waves. When the modulation strength is 
weak, this system is similar to a BEG in PP. When the 
modulation strength is strong enough, even the Bloch 



state at the Brillouin zone center becomes unstable. This 
means that the BEG loses it superfluidity completely. 
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